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Abstract 

Einstein's field equations with cosmological constant are anal- 
ysed for a static, spherically symmetric perfect fluid having constant 
density. Five new global solutions are described. 

One of these solutions has the Nariai solution joined on as an ex- 
terior field. Another solution describes a decreasing pressure model 
with exterior Schwarzschild-de Sitter spacetime having decreasing 
group orbits at the boundary. Two further types generalise the Ein- 
stein static universe. 

The other new solution is unphysical, it is an increasing pressure 
model with a geometric singularity. 
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1 Introduction 



Cosmological observations jSJ ^3 gi ye strong indications for the presence 
of a positive cosmological constant, on the other hand the low energy limit 
of supcrsymmctry theories requires a negative sign [21) . 

Therefore it is interesting to analyse solutions to the field equations 
with cosmological constant representing for example relativistic stars. An 
overview of these solutions is given in table [IJ appendix lAl 

With vanishing cosmological constant the first static and spherically 
symmetric perfect fluid solution with constant density was already found 
by Schwarzschild in 1916 [17) . In spherical symmetry Tolman |2U[ and Op- 
penheimer and Volkoff ^1] reduced the field equations to the well known 
TOV equation. For simple equations of state Tolman integrated the TOV 
equation and discussed solutions. Although he included the cosmological 
constant in his calculations he stated that it is too small to produce ef- 
fects. The geometrical properties of constant density were analysed by 
Stephani [Tg] . 

With non-vanishing cosmological constant A ^ constant density so- 
lutions were first analysed by Weyl [55]. The different possible spatial 
geometries were already pointed out and the possible coordinate singular- 
ity was mentioned. More than 80 years later Stuchlik ^1 analysed these 
solutions again. He integrated the TOV-A equation for possible values of 
the cosmological constant up to the limit A < -iirpo, where po denotes this 
constant density. In these cases constant density solutions describe stel- 
lar models. In the present paper it is shown that a coordinate singularity 
occurs if A > 47rp , already mentioned by Weyl [22) . 

If the cosmological constant equals this upper bound the pressure van- 
ishes at the mentioned coordinate singularity. In this case one has to use 
the Nariai solution [101 II 1| to get the metric C 1 at the boundary. For 
larger cosmological constant the pressure will vanish after the coordinate 
singularity. The volume of group orbits is decreasing and there one has 
to join the Schwarzschild-de Sitter solution containing the r = singular- 
ity. Increasing the cosmological constant further leads to generalisations 
of the Einstein static universe. These solutions have two regular centres 
with monotonically decreasing or increasing pressure from the first to the 
second centre. Certainly the Einstein cosmos itself is a solution. 

Solutions with A > 47T/9 have not been analysed in literature so far. 
Thus the described constant density solutions obeying this bound are con- 
sidered to be new, except for the Einstein static universe. 



2 



2 Cosmological TOV equation 



Consider the following line clement 

ds 2 = -le^Mdt 2 + e a «dr 2 + rW, (2.1) 
A 

where G = 1 and k = 8ttX 2 will be used. A = 1 corresponds to Ein- 
stein's theory of gravitation in geometriscd units, roughly speaking A — > 
corresponds to the Newtonian theory [31EJE]- 

The field equations for a perfect fluid are three independent equations, 
which imply energy-momentum conservation. One may either use the three 
independent field equations or one uses two field equations and the energy- 
momentum conservation equation. For the purpose of this paper it is more 
convenient to do the second. Hence consider the first two field equations 
and the conservation equation 

_L e Mr) jL ( r _ re -a(r)\ _ Ae Mr) = g^^M ( 2 .2) 

Xr 2 dr V / 

-i (l + rXv'(r) - e a(r >) + XAe a ^ = 8ttA 2 P(r)e a{r) , (2.3) 

-^l(\P(r)+p) = pi(r). (2.4) 

There are four unknown functions. From a physical point of view there 
are two possibilities. Either a matter distribution p = p(r) or an equation 
of state p = p(P) is prescribed. The most physical case is to prescribe 
an equation of state. Instead, in the third section solutions with constant 
density are investigated. 

Equation l|2.2() can easily be integrated. Set the constant of integration 
equal to zero because of regularity at the centre. Introduce the definitions 
of mass up to r and mean density up to r by 

m(r) =471-/ s 2 p(s)ds, w(r) = — ^r 1 , (2.5) 



respectively. Then the first field equation 12.21) yields to 

e- Q M = l-A2 W (r)r 2 -A-r 2 . (2.6) 

3 

Therefore one may write metric (|2.1ll as 

ds 2 = -\e x ^dt 2 + — 7 —+r 2 dn 2 . (2.7) 
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The function v'(r) can be eliminated from (|2.3[l and (|2.4[) and yields to the 
TOV-A equation [H EH 

p,, , (A4^P(r)+ W (r)-A)( A P(r)+p(r)) 

p W = ~ r 1 — xo n 2 — 7x1 > ( 2 - 8 ) 

1 — A2u>(r)H — A^r 2 

which for A = gives the Tolman-Oppenhcimer-Volkoff equation, without 
cosmological term, short TOV equation |14l I2U) . 

If an equation of state p = p(P) is given [T^], the conservation equa- 
tion l|2.4(l can be integrated to give 

r p ( r ) 2dP 
„( r ) = _ / (2.9) 

where P c denotes the central pressure. Using the definition of m(r) then 
l|2.5|) and l|2.8|) form an intcgro-diffcrential system for p{r) and P{r). Dif- 
ferentiating (|2.5|l with respect to r implies 

W 'M = - Unp{P(r)) - 3w(r)) . (2.10) 
r 

Therefore given p = p(P) equations i|2.8J) and (|2.10J) are forming a system 
of differential equations in P(r) and w(r). In |16j existence and uniqueness 
for a given equation of state for this system was shown. The solution for 
w(r) together with H2.6(l gives the function a(r), whereas the solution of 
P(r) in (|2.9jl defines v{r). 



3 Solutions with constant density 

For practical reasons the notation is changed to geometrised units where 
c 2 = 1/A = 1. Assume a positive constant density distribution p = po = 
const. Then w = ^-po gives l|2.8|l in the form 

_ _ r (4,P(r) + ^-J)(P( r)+w ,) | 
1 — kr l 



where k is given by 



8tt A 



In the following all solutions to the above differential equation are derived, 
see QjmHn]. 
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The central pressure P c = P(r = 0) is always assumed to be positive 
and finite. Using that the density is constant in l|2.9|l metric <|2.7|) reads 

ds ' = - ( Sft P ° ) 2 dt " + T~T^2 + ( 3 - 3 ) 
\P(r)+poJ 1-kr 2 

For non-vanishing k a new coordinate a will be introduced 

k < 0, r = ~p= sinha, (3.4) 

k > 0, r=—= since. (3-5) 

The metric is well defined for radial coordinates r <E [0,f) if A > — 87rpo, 
where f denotes the zero of g rr . If A < — 87rpo the metric is well defined 
for all r. 

Solutions of differential equation (|3.1f) arc uniquely determined by the 
three parameters p , P c and A. Therefore one has a 3-parameter family of 
solutions. 



3.1 Stellar models with spatially hyperbolic geometry 

A < -87T/9o 

If A < — 87rpo then k < and the spatial geometry of (|3.3|l is hyperbolic and 
the differential equation does not have a singularity. The volume of group 
orbits has no extrema, thus metric (|3.3|l has no coordinate singularities and 
is well defined for all r. 

Integration of l|3.1(l yields to 



1 - -r^- - Ccosha 

P(a)=po + A : Pa \ = ■ (3.6) 

6 cosh a — 3 

C is the constant of integration, evaluated by defining P(a = 0) = P c to 
be the central pressure. It is found that 

Pc + pa 

The function P(a) is well defined and monotonically decreasing for all a. 

The pressure converges to — po as a or as the radius r tends to infinity. 
Thus there always exists an ctb or R where the pressure vanishes. Vanishing 
of the numerator in l|3.6() yields to 

cosha 6 = ^ (1- ] , (3.8) 
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the corresponding radius R is given by 



If A x 2 



R 'w^x\\c-A 1 -^) (3 ' 9) 

The new coordinate a simplifies expressions. Therefore the radial coor- 
dinate will only be used to give a physical picture or if it simplifies the 
calculations. 

From equation l|3.8|) one can deduce the inverse function. Then one 
has the central pressure as a function of a\,-, Pc = Pc{&b)- To do this the 
explicit expression for C from (|3.7|l is needed. One finds 

(l - -r>-) (cosh a fc - 1) 

p c =po \ 47r ;°v -• ( 3 - 10 ) 

[} - J -3cosha 6 

The central pressure given by l|3.10|) should be finite. Therefore one 
obtains an analogue of the Buchdahl inequality. It reads 

cosha 6 < \ ( 1 - ) . (3.11) 

3 V inpoj 

Thus there exists an upper bound for ctb for given A and pQ. Use that 
sinh(arccosh(a)) = V 'a 2 — 1, then the corresponding radius R is given by 

i(4 ^_ 

471-po 

Since the cosmological constant is negative the right-hand side of l|3.12(l is 
well defined. Using the definition of mass M = (4ir/3)poR 3 one can rewrite 
H3.12|) in terms of M, R and A. This leads to 



3M<^R + R^-±R2. (3.13) 

This is the wanted analogue of the Buchdahl inequality l|3.27f) . 

At r = i?, where the pressure vanishes, the Schwarzschild-anti-de Sitter 
solution is joined. 



Joining interior and exterior solution 

At the P — surface the Schwarzschild-anti-de Sitter solution is joined by 
using Gauss coordinates to the r = const, hypcrsurfaces and defining the 
mass by M — (47r/3)po-R 3 - If necessary one can rescale the time to get 
the <7tt-component continuous. With this rescaling and the use of Gauss 
coordinates one can get metric C 1 at the boundary. 

Since the density is not continuous at the boundary the Ricci tensor is 
not, either. Therefore the metric is at most C . This cannot be improved. 
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3.2 Stellar models with spatially Euclidean geometry 

A = -87T/9o 

Assume that cosmological constant and constant density are chosen such 
that 87rp = —A. Then k = and the denominator of (|3.1[) becomes one 
and the differential equation simplifies to 

^ = -47TT (P + po)\ (3.14) 

and the t = const, hypersurfaces of l|3.3|) arc purely Euclidean. As in the 
former case metric l|3.3|) is well defined for all r. Integration yields to 

P(r) = ~ * i -Po, (3-15) 

where the constant of integration is fixed by P(r = 0) = P c . 

The denominator of the pressure distribution cannot vanish because 
central pressure and density are assumed to be positive. Therefore (|3.15|l 
has no singularities. As the radius tends to infinity the fraction tends to 
zero and thus the pressure converges to — po- 

This implies that there always exits a radius R where P(r = R) = 0. 
Therefore all solutions to (|3.14|l in (|3.3|l describe stellar objects. Their 
radius is given by 

^ = ^(--^—)- (3-16) 
2tt \p Pc + PaJ V ; 



One finds that the radius R is bounded by 1/ ^/2irpQ. Inserting the 
definition of the density yields to 

9 1 1 AwR 3 

R 2 < = , 3.17) 

2np Q 2tt 3M ' V ' 

which implies 

3M < 2R. (3.18) 

Thus 1|3.18|) is the analogous Buchdahl inequality to equation l|3.27|) and 
equals (jXT3|) with A = -87rp = -6M/R 3 . 

At r = i?, where the pressure vanishes, the Schwarzschild-anti-dc Sitter 
solution is joined uniquely as described. Since the density is not continuous 
at the boundary the Ricci tensor is not, cither. Thus the metric is again 
at most C 1 . 



7 



3.3 Stellar models with spatially spherical geometry 

A > -87T/90 

If A > — 87rpo then k is positive. Equation 1(3. 1|) has a singularity at 

r = r=-L. (3.19) 

In a the singularity f corresponds to a(f) = 7r/2, the equator of a 3-sphcrc. 
The spatial part of metric 1(3. 3|1 now describes a part of the 3-sphere of 
radius 1/y/k. The metric is well defined for radii less than r. 
Integration of 1)3.1(1 gives 

(t^-I) +Ccosa 

P(a)=Po^—l , (3.20) 

3 — (J cos a 

where the constant of integration C is given by (|3.7|l . One finds 

^/^fG^- 1 )' (3 - 2i) 

which is less than zero if A < 47rpo, which is the restriction in |19| and the 
output of [121 113| . Since this is the considered case the singularity of the 
pressure gradient is not important yet. 

Since P(tt/2) < there exists an at such that P(at) = 0. at can be 
derived from equation l|3.20|) and one obtains the analogue of l|3.9|) 

00805 = 1(1-^). (3.22) 

It remains to derive the analogue of the Buchdahl inequality for this 
case. One uses (13.2211 to find P c = P c (at) 

1 - t^-) (1 - cosa fc ) 

Pc = Po " " P0) f (3-23) 

3 cos at — ( 1 — i^— 

which is similar to 1)3.9(1 . Finiteness of the central pressure in ((3.23(1 gives 
CO sa fe >i(l-^_). (3.24) 



With sin(arccos(a)) = vl — a 2 , this yields to the wanted analogue 



1/4. A 



2 . 3 V 4irp 



R Z < ^ ^, (3-25) 
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which equals (|3.12[1 . The greater sign reversed because (arccosa) is de- 
creasing. Thus it can also be rewritten to give (|3.13[l . At a = or r = R 
the Schwarzschild-anti-de Sitter solution is joined and the metric is C 1 . 
Without further assumptions this cannot be improved. 



3.4 Solutions with vanishing cosmological constant 

A = 

Assume a vanishing cosmological constant. Then one can use all equations 
of the former case with A = 0. Only one of these relations will be shown, 
namely the Buchdahl inequality It is derived from H3.25|> 

R 2 < (3.26) 

using M — (4tt/3)poR 3 leads to 

M < ^-R, 3M < -R. (3.27) 
9 3 

3.5 Stellar models with spatially spherical geometry 

< A < Anpo 

Integration of the TOV-A equation gives H3.20|) . Because of l|3.21|l the 
boundary P(R) = exists. As in the former sections one finds (|3.12|l and 
written in terms of M, R and A again gives 



3M< 1r + R^-±R2. (3.28) 

Since the cosmological constant is positive the square root term is well 
defined if 

R<f^- (3.29) 

In this section < A < 47rpo is assumed. Using the definition of mass this 
can be rewritten to give 



i? 2 A < 4^ i? 2 = ^ < | + - ±R\ (3.30) 

where the right-hand side of (|3.28|l was used. By taking the 2/3 on the 
other side and squaring the equation this reduces to the simple inequality 

i? 2 A<l, R<-=, (3.31) 
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which means that the boundary of the stellar object is located before the 
cosmological event horizon. 

It is remarkable that bounds on mass and radius similar to (|3.28[) 
and (|3.31[) were derived in |12l I13| by considering the virial theorem in the 
presence of a positive cosmological constant. Furthermore the inequality 
A < 4irp is derived, which in the context of the this paper may be seen as an 
evidence that the exterior spacetime of stellar objects is the Schwarzschild- 
dc Sitter spacetime with increasing group orbits at the boundary. 

As before the Schwarzschild-de Sitter solution is joined C 1 . Since the 
density is not continuous at the boundary this cannot be improved. 

One can construct solutions without singularities. The group orbits are 
increasing up to R. The boundary of the stellar object r = R can be put 
in region / of Pcnrose-Carter figure ^ where the time-like killing vector 
is future directed. This leads to figure the spacetime still contains an 
infinite sequence of singularities r = and space-like infinities r = oo. It 
is possible to put a second object with boundary r = R in region IV of 
Penrose- Carter figure 02 where the time- like killing vector is past directed. 
This leads to figure 0] and shows that there are no singularities in the 
spacetime. This spacetime is not globally static because of the remaining 
dynamical parts of the the Schwarzschild-de Sitter spacetime, regions IIq 
and I lie- Penrose-Carter diagrams are in the appendix 151 

A black hole alternative solution without singularities and without hori- 
zons was described by [U]. 

3.6 Solutions with exterior Nariai metric 

A = 4np 

In this special case integration of (|3 . 1|> gives 

^/ \ Ccosa ,„ „„. 

P « = Po^^ . 3.32 

3 — G cos a 

The pressure ()3.32ll vanishes at a = at, = ir/2, the coordinate singularity 
in the radial coordinate r, the interior metric reads 

One would like to join the Schwarzschild-de Sitter metric, but this does 
not work, since the group orbits of the Schwarzschild-dc Sitter metric are 
increasing whereas the group orbits of the interior metric have constant 
volume. But there is one other spherically symmetric vacuum solution to 
the Einstein field equations with cosmological constant, the Nariai solution 
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[TOl EH- Its metric is 



d « 2 = A 



(A cos log r + Bsinlogr) 2 dt 2 + \dr 2 + dO 



where A and B are arbitrary constants. With r = e a this becomes 
ds 2 = - [-(Acosa + Bsina) 2 dt 2 + da 2 + dft 2 ] . 



(3.34) 



(3.35) 



Metrics (|3.33|l and (|3.35|l can be joined by fixing the constants A and B. 
With 



Po 



b= Pc + Po Va, 

Po 



(3.36) 



the metric is C 1 at a = ir/2. As before, since the density is not continuous, 
the metric is at most C 1 . 



3.7 Solutions which have decreasing group orbits at 
the boundary 

47rp < A < Ao 

Integration of the TOV-A equation again yields (|3.20J) . Assume that the 
pressure vanishes before the second centre a = tt of the 3-sphere is reached. 
The condition P(a = tt) < leads to an upper bound for the cosmological 
constant. This bound is given by 

mm- 

Then iirpo < A < A implies the following: 

The pressure is decreasing near the centre and vanishes for some ctb, 
where ir/2 < < tt. Equation 1)3. 5|l together with the metric (|3.3|) implies 
that the volume of the group orbits is decreasing if a > ir/2. 

At ah one uniquely joins the Schwarzschild-de Sitter solution by M = 
(4tt/3)poR 3 . With Gauss coordinates relative to the P(«b) = hypersur- 
face the metric will be C 1 . But there is a crucial difference to the former 
case with exterior Schwarzschild-de Sitter solution. Because of the de- 
creasing group orbits at the boundary there is still the singularity r = 
in the vacuum spacctimc. Penrosc-Carter diagram [S] shows this interesting 
solution. 
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3.8 Decreasing solutions with two regular centres 

A < A < A E 

As before the pressure is given by (|3.2U[) . Assume that the pressure is 
decreasing near the first centre a = 0. This gives an upper bound of the 
cosmological constant 

A E :=47rp (3P c /po + l), (3.38) 

where A^ is the cosmological constant of the Einstein static universe. These 
possible values A < A < A^ imply: 

The pressure is decreasing is near the first centre a = but remains 
positive for all a because A > Aq. Therefore there exists a second centre 
at a = 7r. At the second centre of the 3-sphere the pressure (|3.2U[) becomes 

P (<* = «)=P°^ f Y T ^ ■ (3-39) 

It only vanishes if A = A . The solution is inextendible. The second centre 
is also regular. This is easily shown with Gauss coordinates. 

Solutions of this kind are generalisations of the Einstein static universe. 
These 3-spheres have a homogeneous density but do not have constant 
pressure. They have a given central pressure P c at the first regular centre 
which decreases monotonically towards the second regular centre. Gener- 
alisations of the Einstein static universe were also found earlier in Rcf. 7 
by Ibrahim and Nutku. 1 

3.9 The Einstein static universe 

A = Ae 

Assume a constant pressure function. Then the pressure gradient van- 
ishes and therefore the right-hand side of (|3.1|l has to vanish. This gives 
an unique relation between density, central pressure and cosmological con- 
stant. One obtains (|3.38(l 

A = A E = 4tt (3P e + Pa) , (3.40) 

where P E = P c was used to emphasise that the given central pressure 
corresponds to the Einstein static universe jH] and is homogeneous. Its 
metric is of no further interest. 

For a given density pa there exists for every choice of a central pressure 
P c a unique cosmological constant given by (|3.38J) such that an Einstein 
static universe is the solution to H3.1|) . 

I would like to thank Aysel Karafistan for bringing this reference to my attention. 
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3.10 Increasing solutions with two regular centres 

A E < A < A s 

Assume that the pressure (13.2011 is finite at the second centre. This leads 
to an upper bound of the cosmological constant defined by 

As :=4 7 r / 9 (6P c / / 9o+4). (3.41) 

The possible values of the cosmological constant imply: 

The pressure P(a) is increasing near the first regular centre. It increases 
monotonically up to a — tt, where one has a second regular centre. This 
situation is similar to the case where Ao < A < Ae- These solutions are 
also describing generalisations of the Einstein static universe. 

The generalisations are symmetric with respect to the Einstein static 
universe. By symmetric one means the following. Instead of writing the 
pressure as a function of a depending on the given values po, Pc and A one 
can eliminate the cosmological constant with the pressure at the second 
centre P(a = tt), given by H3.39(l . Then one easily finds that the pressure is 
symmetric to a = tt/2 if both central pressures are exchanged and therefore 
this is the converse situation to the case where Ao < A < Ae- If the central 
pressures are equal the dependence on a vanishes and one is left with the 
Einstein static universe. 

3.11 Solutions with geometric singularity 

A> A s 

In this case it is assumed that A exceeds the upper limit Ag. Then (|3.2U|I 
implies that the pressure is increasing near the centre and diverges before 
a = tt is reached. Thus these solutions have a geometric singularity with 
unphysical properties, as can be seen by considering the square of the 
Riemann tensor. Therefore they are of no further interest. 
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A Overview of constant density solutions 



cosmological 
constant 


spatial 
geometry 


short description of the solution 


A < -871-po 


hypcrboloid 


stellar model with exterior 
Schwarzschild-anti-de Sitter solution 


A = -871-po 


Euclidean 


stellar model with exterior 
Schwarzschild-anti-de Sitter solution 


-8irp < A < 


3-sphere 


stellar model with exterior 
Schwarzschild-anti-de Sitter solution 


A = 


3-sphere 


stellar model with exterior Schwarzschild 
solution 


< A < Anpa 


3-sphcrc 


stellar model with exterior 
Schwarzschild-dc Sitter solution 


A = AnpQ 


3-sphere 


stellar model with exterior Nariai solu- 
tion 


Airpo < A < A 


3-sphere 


decreasing pressure model with exte- 
rior Schwarzschild-de Sitter solution; the 
group orbits are decreasing at the bound- 
ary 


A < A < A E 


3-sphere 


solution with two centres, pressure de- 
creasing near the first; generalisation of 
the Einstein static universe 


A = Ae 


3-sphere 


Einstein static universe 


A e < A < A s 


3-sphere 


solution with two centres, solution in- 
creasing near the first; generalisation of 
the Einstein static universe 


A>A S 


3-sphere 


increasing pressure solution with geomet- 
ric singularity 



Table I: Overview of constant density solutions 



14 



B Penrose- Carter diagrams 




Figure 1: Penrose-Carter diagram for Schwarzschild-de Sitter space. The 
Killing vector d/dt is time-like and future-directed in regions / and time- 
like and past-directed in regions IV. In the others regions it is space-like. 
The surfaces r = r + and r = r ++ are black-hole and cosmological event 
horizons, respectively. and are the space-like infinities. 




Figure 2: Penrose-Carter diagram for Schwarzschild-anti-de Sitter space. 
The surface r = r + is the black-hole event horizons. and ^~ are the 
time-like infinities. 
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r = oo, r = 




r = oo, r = 



Figure 3: Penrose-Carter diagram with one stellar object having a radius 
R which lies between the two horizons. The group orbits are increasing 
at the boundary. The shaded region is the matter solution with regular 
centre. There is still an infinite sequence of singularities r = and space- 
like infinities r = oo. 



r = oo, 




r = oo, J? 



Figure 4: Penrose-Carter diagram with two stellar objects having radii R 
which lie between the two horizons. Since the group orbits are increasing 
up to R the vacuum part contains the cosmological event horizon r ++ . This 
solution with two objects has no singularities. Because of regions lie an d 
IIIc this spacetime is not globally static. 
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Figure 5: Pcnrosc-Carter diagram with two stellar objects having radii 
R which lie between the two horizons. The group orbits of the interior 
solutions are decreasing where the vacuum solution is joined. Thus the 
r = singularity of the vacuum part is present. 
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